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We compute the gravitational waveform produced by cosmic superstring re- 
connections. This is done by first constructing the superstring reconnection 
trajectory, which closely resembles that of classical, instantaneous reconnection 
but with the singularities smoothed out due to the string path integral. We then 
evaluate the graviton vertex operator in this background to obtain the burst 
amplitude. The result is compared to the detection threshold for current and 
future gravitational wave detectors, finding that neither bursts nor the stochastic 
background would be detectable by Advanced LIGO. This disappointing but an- 
ticipated conclusion holds even for the most optimistic values of the reconnection 
probability and loop sizes. 
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1 Introduction 



One of the most exciting products of the recent synthesis of superstring theory and 
cosmology has been that of cosmic superstrings [1] . Like their classical cousins [2] , cos- 
mic superstrings were briefly considered then discarded [3] , but for theoretical reasons 
rather than observational ones. Subsequent (and largely non-perturbative) analysis [4] 
has shown that cosmic superstrings have been found to naturally arise in models of 
brane inflation [5] , a string theory realization of inflationary cosmology theory. 

Cosmic strings and superstrings can produce a variety of astrophysical signatures 
including gravitational waves [6-10], ultra high energy cosmic rays [11], and gamma 
ray bursts [12]. Very recent work has revealed a number of exciting new possibilities 
such as radio bursts from strings [13], effects on the cosmic 21 cm power spectrum [14], 
magnetogenesis [15], effects on the CMB at small angular scales [16,17], CMB po- 
larization [18], microlensing from strings [19,20], strong lensing [21,22], and weak 
lensing [23]. 

Prior work on gravitational radiation has focused on two processes: cusps (whereby 
a segment of the string momentarily moves at the speed of light) and kinks (formed 
after two cosmic strings collide and reconnect) [6-10]. Here we study a third source: 
the radiation emitted from the reconnection process itself. This is possible for cosmic 
superstrings because string theory allows us to explicitly construct the reconnection 
process and compute detailed interaction properties. We will do this for the bosonic 
string but the presence of fermions in the superstring should not change the conclusions. 

In the first half of this article we compute the gravitational waveform resulting 
from fundamental cosmic superstring reconnection. In the second half we consider the 
likelihood of detection of this signal with current and future experiments. 

2 Cosmic Superstring Reconnection 
2.1 Basic Process 

The theory of cosmic superstring reconnection was developed in [24] [25] [26]. Consider 
two long, straight wound bosonic strings on a 2D torus of size L and skew angle 9 as 
illustrated in Figure 1. In terms of the string tension (27rQ;')~^ the momenta are taken 
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Figure 1: We model cosmic superstrings as straight wound modes on a large torus, 
which will then interact to form a kinked configuration. 
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(1,0,0,0,0), Li = L(0,l,0,0,O), (2.1) 

-^/^{l, 0, 0, V, 0), L2 = L(0, cos 9, sin 9, 0, 0). 



These satisfy the tachyonic mass-shell conditions 
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The relevant vertex operators are: {i — 1,2) 



VT{z,z;pi) 



(2.2) 



2^:^/V 

where the volume V — V^L'^ sin 9 is the product of the the transverse volume and the 
2D torus (methods to calculate V±_ can be found in [27]). Here Xl{z),Xr{z) refer to 
the (anti)holomorphic components of X{z, z), so that X{z, z) — Xl{z) + Xr{z). 

These will scatter into some kinked configuration and the amplitude for all such 
processes can be summed. In the large-winding limit the probability of reconnection 
is found to be 
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It is also likely that there will be radiation emitted during this reconnection. This 
can be included by using the reconnection process as a classical background trajectory 
Xci{z, z) upon which the radiation vertex operator must be integrated over [26]. Then 
the probability to emit a radiated state of definite momentum k will be the probability 
of reconnection times a factor depending on the radiated particle under consideration, 



One might also consider the possibility that the strings will pass through each other 
without reconnecting but still emit radiation in the process. This is a valid possibility 
but one in which the amplitude is suppressed by a factor of Qs due to three final 
strings instead of two (recall that a reconnected kinked string is considered a single 
string, whereas this process will produce two unconnected strings with accompanying 
coupfing factor). So to leading order this process can be neglected. 

2.2 Computing the Background Trajectory 

The reconnection trajectory was computed in [26] by considering the vertex operator 
of the kinked string. This is done by examining the operator product expansion of the 
straight strings. 



The Taylor expansion of the exponential shows the vertex operators of the infinite 
tower of the produced states, which will appear kinked due to their large oscillator 
excitation number N: 



Since pm ■p2R = PiL ■P2L the result will be identical for the right-moving oscillators and 
so N — N. Now consider the state corresponding to the sum of these vertex operators 




(2.4) 
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Ikinks) = -i^ /rfede |e|-W)\/T(e,e;Pi)^T(0;p2) 



dede 



\Pl+P2)- 



(27r2)2 

The expectation value of the position operator in the reconnection amphtude is then 
easily evaluated: 

(Vkinks(oo; pi,p2)X{z, z)Vt{1;pi)Vt{0; P2)) 



Xr] (Z, Z) 



(Vkinks(oo; Pi,P2)Vt{1;Pi)Vt{0; P2)) 
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While this is a technically correct answer, it is unsatisfactory in that it appears to treat 
string 1 different from string 2 even though of course these are unphysical labels. The 
reason is that the placement of vertex operators has treated them in a non-symmetric 
fashion, expanding Vrie; pi) around the location of Vt(0; ^2)- To remedy this we simply 
redo the previous calculation but expand both vertex operators around the midpoint 
of their separation. 



Ikinks) = ^^^^dede |e|-W)v^(e/2;p,)Vy(-e/2;p2) 



(27ri) 



(2.5) 



dede |e| -^(^+^)eV ^">l Pli-«-n(e/2)Vn+P2L-a-n(-e/2)Vn+{i^-R) _^ _ 



A good check that this is a physically equivalent solution (i.e. obtainable via a confor- 
mal transformation z'{z)) is that for the massless kink state {N = 1) the polarization 
vector piL,B. has merely been replaced with {pil,r — P2L,_r)/2, which is a difference 
proportional to the wavevector Pil,r + P2L,r and thus pure gauge. Now combining 
this new kink vertex operator with symmetric relocations 2 = ±| for the incoming 
straight-string vertex operators, the transformed trajectory is 

(^kinks(oo;pi,p2)^(2,^)V'r(|;pi)VT(-^;P2)) 
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As hoped, this is symmetric under pn ^ p2L by exchanging z ^ —z (up to sub- 
tleties involving branch cuts). This simplifies in the cosmic string limit (A^ oo) by 
performing the summation, 
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We will neglect the second term since this is a good approximation for most of the 
parameter space, and the trajectory simplifies tremendously: 
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+ {L ^ R,z ^ z). (2.6) 



However, the trajectory (2.6) as written is complex, due to the analytic continuation 
required by string theory; we now need to transform this back into a trajectory in real 
space. This is easily done by recalling that vertex operators are inserted into the path 
integral as source terms via 
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This produces a trajectory 
X(z,z) 



^ j |(9X|2 - ip^XLizo) - ipRXji{zo) 



-iyPi \n{z 

I 

.Oi 



L 



—i — p In I z — 1 —i — In 



zo 



AtT Z — Zq 

We see the winding term is real, and that moving around the vertex operator produces 
X X + L as desired. The momenta term, however, is purely imaginary; this can 
be remedied by rotating p ip, allowing us to interpret the (logarithmic) distance 
away from the vertex operator as displacement in time. Thus the trajectory for the 
reconnection process in real target space is 
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(2.7) 
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This reconnection process is shown in Figure 2, from the worldsheet and spacetime 
viewpoints. The infinite past (containing the two incoming asymptotic states of simple 



6 



wound strings) are located at z = the infinite future (containing tlie single outgoing 
asymptotic state of a sharply kinked reconnected string) is located at 2; = ±2, and the 
entire unit circle |z| = 1 maps to time X° = 0. As expected, string theory has smoothed 
out what is classically a very violent and singular process. 

2.3 Gravitational Radiation 

Now let us compute the radiation emitted from such a reconnection process. The 
simplest and most physically relevant vertex operator is that for a graviton, 

^grav = e^,9X£9X]^e-*'=-^^/=^-*'=-^«/^ (2.8) 

where we have factored the left- and right-moving components. Our task now is to 
evaluate the amphtude 

Aik) = ^ Jdh Vg,^,[X,i{z, z)] (2.9) 



given the trajectory (2.7). As explained in more detail in section 3.1, this problem is 



tantamount to determining the cosmic string's energy-momentum tensor 

T^^'ik) = — [ (fz dXldX^'^e-'^-^^/^-'^-^^/^. (2.10) 
2Txa' J 

Damour and Vilenkin investigated this quantity for classical cosmic strings [6], finding 
that it is largest when both the left- and right-moving components of the exponential 
have either a saddle point or discontinuity at the same location. The signal strength 
can then be conveniently ranked according to the phase behavior: 

1. Cusp: double saddle-point. 

2. Kink: one saddle-point, one discontinuity. 

3. Reconnection: double discontinuity. 

Since the contribution of radiation emitted from reconnection would appear to be a 
sw&-sub-leading correction, it is not surprising that henceforth no effort has been made 
to explicitly calculate this. 

The reason we believe this is now a worthwhile endeavor is because in superstring 



theory the probability of reconnection (2.3) can be greatly suppressed, P ~ 0.01. As 
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Figure 2: String reconnection from the worldsheet and spacetime perspectives. The 
upper panels show the equal-time contours before, during, and just after the recon- 
nection process. The lower panels show the corresponding spacetime trajectories. We 
have labeled the portion of the contours in red (blue) corresponding to the asymptotic 
first (second) strings, but of course during reconnection the distinction is irrelevant 
since it is now one large string. g 



discussed in more detail in the next section, this then produces a larger number of 
strings ~ 1/P, the number of reconnection attempts now scaling as ~ 1/-P^, and the 
number of successful reconnections scaling as ~ Thus, as is the case for bursts 

from cusps, the lower the probability of reconnection, the larger the gravitational wave 
signal. 



Let us make a naive attempt to evaluate the amplitude (2.10). In examining X^i in 



(2.7), there do not appear to be any kink-like discontinuities, but there do appear to be 
cusps as indicated by the presence of extrema. Focusing now on just the holomorphic 
part of Xci, the cusp condition is 

Rewriting this as a quadratic equation, one sees there are two solutions, 

5{l + rL)±^/i9 + rL){l + 9rL) k ■ p^L 

4(rL - 1) k- piL 

One solution will be ~ and the other ^ oo. In the large- winding limit L oo, 
away from these extrema the phase will oscillate wildly and so we attempt to utilize a 
saddle-point approximation. 



T^^ik) 



^k ■ d^X^iz^ ^ ^d^XLiz^) 



x[L ^ R,z ^ z) (naive) 



This could not be possibly be correct, as this T L/k spectrum indicative of si- 
multaneous cusps does not even remotely reproduce the T ~ Ijk"^ spectrum due to 
simultaneous kinks, which we must obtain in the classical limit. Where did we go 
wrong? The answer lies in the fact that our saddle-point evaluation was done at the 
edge of moduli space 2;+ ~ oo, where such an approximation breaks down [28]. 

To evaluate the amplitude correctly, we follow Kawai-Lewellen-Tye [29] and write 
the exponential part of the amplitude explicitly: 



A{k) ~ / (i'z 



Is \k-p^L / ^ I 1 \ ^k-p2L / ^ _ 1 \ / ^ I 1 \ ^k-p2R 

2 \ I ^ 2 \ I 2 \ ~^ 2 



- Vl + -/ \1 - - J V1 + 

2/ \-^^2/ \ 2/ \-^^2 



Denoting z = x + iy, the integrand is an analytic function of y with eight branch 
points: ±i (a; ± ^) , ±i (a; ± 2). The contour over y stretches from —00 to 00 along the 
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Figure 3: The contours for z^, initially lie on the real axis (above) and are chosen to 
circumnavigate the poles so as to produce the correct phases in the amplitude. These 
can then be deformed into something more convenient for explicit evaluation, such as 
the complex unit circle with negligible contribution from infinity (below). 



real axis, so we can then rotate it counter-clockwise to now run from —oo to oo along 
the imaginary axis. If we then define the independent variables 



zl = x + iy, 



ZR = x-iy 



we can rewrite the integral in terms of these variables, paying special attention to the 
phases. 
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An examination of these phases produces the following summary: 
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Figure 4: The reconnection geometry for classical cosmic strings (black) versus cosmic 
superstrings, both exactly (red) and in the large-L approximation (blue). 
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These phases can be produced by choosing the contours shown in Figure 3. Once 
these choices have been made, of course, we are free to deform them into something 
more convenient for analysis. We choose to deform the contour into the complex unit 
circle, which we recall corresponds to zero time as shown in Figure 2. Note that the 
trajectory at this moment resembles a smoothed-out form of the classical reconnection 
picture. The idea of transforming an amplitude at a given point on the string over all 
time into over the entire string at a single moment in time is not without precedent in 
superstring theory [30]. 

Now that we have completely factored the left- and right moving components to 



(2.10), we are only concerned with the evaluation of the integral over the angle 
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arg zl, 



(2.11) 



Defining the scale-invariant quantities = Xi/L^ipi = Lcf)^, we can expand (2.11) 



in powers of L near the points 2;/, = ±1 at which the trajectory is very nearly linear: 
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(2.12) 
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Thus in accordance with the classical expectation, we have obtained a spectrum which 
is scale-invariant in the large- winding limit, and which has virtually identical depen- 
dence upon the relative angle and velocity of the strings. In the last line we have 
written it to most resemble the classical answer plus corrections; the terms propor- 
tional to ^ are the corrections from the string path integral. Of course the procedure 
is then repeated for the right-moving side. 

We must also average over all geometric possibilities for which reconnection takes 
place using the probability measure 



dP 



dv^ 
J -I 



1 ^27r 

(i(cos 9) 



P{v,9) 



(2.13) 



modulated by some string distribution n{v, 6) and then finally to sum over polariza- 



tions. While (2.13) itself easily evaluates to / dP 



averaging (2.12) must be 



done numerically. 



2.4 Other Cosmic Superstring and Radiation Types 

Though we have considered massless gravitational modes, the same procedure will work 



for any string mode by using the appropriate vertex operator in (2.9 ). However, analysis 
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of the decay modes for highly-excited strings suggests that radiation into massive states 
is suppressed and that only massless quanta are excited [31]. 

In addition to the fundamental cosmic superstrings considered here, there are also 
D- and {p, g)-strings produced in brane inflation which will have significantly larger ten- 
sion are therefore might be capable of emitting higher-energy bursts [32] [33] . Though 
the method of calculation for reconnection interactions is very different [25] [26] [34] 



[35], ideally their waveform could be estimated by adapting the F-string result (2.12) 



and appropriately replacing the inverse tension as a' — > a'/yp^ + q'^/g'^ in the BPS 
limit, depending on the specific brane inflation model. These {p, q) strings also form 
more sophisticated networks [36] and thus would have distinctive burst patterns [37]. 

2.5 Classical Cosmic String Reconnection 

For verification and comparison we now compute the string trajectory of two classical 
Nambu-Goto strings during reconnection. This is done using the geometry in Figure 1, 
taking the point of reconnection to be cr = and the gauge conditions for the spatial 
coordinates X(r, a) to be 

X'2 + X2 = 1, 

whereas the time coordinate is defined to be = 1. Though there exist analytic and 
numerical studies of the reconnection process [38] [39] [40], here we take the simplistic 
choice of instantaneous reconnection. At time r = we can write the spatial derivative 
and velocity for the first string as 

X cr < 

■^^^{cosipx + sinipy) cr > ' 



X'l(0,(7) = 

Xi(0,(T 



a < 

I3z a>0 ' 



and for string 2 as 

X'2(0,a) 



7 ^{cosipx + sintpy) cr < 
X cr > 



X2(0,(t) 



I3z a <0 
a > 
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Now defining cr-t = r ± cr, we can always write X(r, cr) = X+(o"+) + X_(o"_) and 
then these tangent vectors combine into 



5±Xi(a±) 



9±X2(a±) 



±X ±(T± < 

[3z ± 7^-'-(cos'?/'x + sinT/'y) ±(T-|- > ' 

[3 z {cos Tpx + sm'ipy) ia-t < 

±x ±a± > 0. 



(2.14) 



These are then utihzed in the evaluation of the integral (2.11), which is easily done 
near a discontinuity: 



2i 



da± d^X^'e-'^-^^ 
9X^^(0+) 9X^(0-) 



To compare these to the superstring results, recall that in (2.1) we have defined the 
left/right momentum vectors Pil,r,P2L,r, which can then be normalized as (substitut- 
ing + ^ L, > R), 



P'± 



^ 2-Ka' 



The classical result (2.14) can then be directly compared to the superstring result (2.12) 
in terms of the geometry at the moment of reconnection, as shown in Figure 4. 



3 Detectability of bursts and stochastic background 

In a network of cosmic superstrings reconnections can occur in two situations, (1) when 
two long strings meet, and (2) when loops are formed. The rate at which reconnections 
take place in a network of superstrings, and the gravitational wave bursts that result, 
can be readily estimated. 

3.1 Cosmic Superstring Detectability 

If we assume the string density to be inversely proportional to the reconnection prob- 
ability, at cosmic time t the density of long strings in the scaling regime is [2,5,25,38] 
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where /x is the mass per unit length in strings and P is the reconnection probabihty. 
The total energy in long strings is ~ pt^ ^ fit/P and the total length is ^ t/P, so the 
total number of long strings in a Hubble volume at any time is ~ P~^ [41]. Strings move 
relativistically so that in a time t they travel a distance ~ t and each string attempts 
to reconnect with every other string in a Hubble volume. Although the number of 
reconnection attempts in a time t is ~ P~^ the number of successful reconnections is 
~ P~^. So the number of long string reconnections that take place per unit space-time 
volume is 

^R,s-j^- (3.15) 

Reconnections also occur when loops are formed. The scaling network converts a 
length of string ~ t/P into loops in a Hubble time ~ t per Hubble volume (~ t^). If 
loops have a size at at formation, the number of loops formed in a time t is ~ {aP)~^, 
so the number of reconnections per unit space-time volume is 

nn,L - J^. (3.16) 



If a is comparable to the Hubble length then the loop contribution Eq. (3.16) to the 



burst rate is comparable to the string contribution Eq. (3.15), but if the size of loops 
a is determined by gravitational back-reaction then the loop contribution dominates. 
The total number of reconnections per unit spacetime volume is thus, 

nn = nn,s + nn,L ^ ]^ + ^ (3-17) 

if a < 1. 

The gravitational wave strain produced at the reconnection event is of order 

Gfi 

so that for a reconnection occurring at a redshift z, 

H()Gjj, 



h{f,z) 



The amplitude of the waveform is suppressed by an additional factor of (1 + z) due to 
the redshifting of the frequency / (see [6]). The Hubble parameter today is Hq = 73 
km/s/Mpc [42]), and we have written the proper distance as r{z) = HQ^Lpr{z). 
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Unlike cusps, these bursts do not point in any particular direction. However, we do 
not observe bursts at arbitrarily low frequencies, rather, only down to the inverse of 
the scale at which the strings curve (~ at) and the behaviour of the waveform is 
no longer valid. Using Eq. (3.17) we can write the number of bursts originating from 
a volume at redshift z as 



dV{z) " aPffiz) 

where the factor of {1 + z)~^ comes from the relation between the observed burst rate 
and the cosmic time, and the 0-function cuts off events that we would not observe 
below frequency /. We also take the volume element at a redshift z to be dV{z) = 
H^^ipy[z)dz and so the rate of events from redshifts between z and z + dz is 

^ ~ ^° iVv^4|~)V v-(^)e((l + z)faH^'v,{z) - 1). (3.18) 

To compute the observable rate of bursts we use the similar methods to compute 
sensitivity estimates as used in [7]. We write the strain produced by the reconnection 
cusp as 

h{f,z) = A{z)r\ 

with 

A{z) ~ (3.19) 

Using the noise curve Eq. (39) of [7] we find a minimum detectable amplitude Amin = 
10-19 s for Initial LIGO and A^^^ = IQ-^o s for Advanced LIGO. Using Eq. ( |3.19[ ) we 
can find what redshifts these minimum amplitudes correspond to for a given value of 
Gjj. (-2max, say) and integrate Eq. ( |3.18[ ) to find the rate 

R ~ r^'^ dzHo ^^^l] ' m^M^ + z)faH^'^,{z) - 1). 
Jo (^Pftiz) 

The cosmological functions that enter these expressions, (pt{z), (pr{z), and (pv{z) 
depend on cosmological parameters and need to be computed numerically as described 
in Appendix A of [7], where a vanilla A-CDM model was used [42]. However, one 
can obtain fairly good analytic approximations for the cosmological functions using 
variations of the expressions introduced in [6] . In particular the following expressions, 

ift = {i+zr'/\i+z/z,^r'/', 
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with Zeq = 5440, 

iPr = ^(1 + ^/3.5)~"'", 

and 

agree to better than 20% with numerically computed functions at all relevant redshifts. 
We will use these functions to compute the burst rates and stochastic backgrounds. 

Figure [5] shows the number of burst events per year as a function of string tension 
for for Advanced LIGO as a function of string tension for a reconnection probability 
P = 10^^. We have taken the minimum detectable amplitude to be Amin = 10~^° s, 
set the frequency to / = 75 Hz, and taken F = 50 everywhere. The solid, dotted and 
dashed-dotted lines show the number of events for loop sizes given by gravitational 
back-reaction [43-47] a = eTGfi with e = 1, e = 10^^, and e = 10^^^ respectively. 
The dashed line shows the burst rate for large long lived loops (a = 0.1) as suggested 
by recent numerical simulations [48] . Unfortunately, none of these models result in a 
signal detectable by LIGO. 

To compute the stochastic background of gravitational waves generated by all the 
reconnections occurring in the network we can use the results in [9] for cusps. The 
spectrum of gravitational waves is given by 

A-tt"^ r°° rl¥} 
figw(/) = ^f/ dzh\f,z)— (3.20) 

^ ^min 

where Zrmn excludes large bursts that occur infrequently and contaminate the estimate 
of the strain [6,9], and h{f, z) is the strain produced at frequency / by a burst from a 
reconnection at redshift z. 

Figure [6] shows three examples of the value of f^cw at the frequencies of various 
experiments as a function of the string tension for P = 10~^. The solid lines show 
the stochastic background generated at 75 Hz, corresponding to LIGO frequencies, for 
two values of the loop size. The dotted line shows the spectrum at 10~® Hz, the band 
at which pulsar timing experiments are most sensitive. Again, unfortunately, none 
of these models result in a detectable signal. Although not shown here, large loops 
(a = 0.1) do not result in a detectable signal. 
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Burst rate for Advanced LIGO, P=10"^ 

10° I r 




Figure 5: Plot of the number of burst events per year for Advanced LIGO (^min = 
10"^'^ s, / = 75 Hz) as a function of string tension for a reconnection probability 
P — 10~^. The solid, dotted and dashed-dotted lines show the number of events for 
loop sizes given by gravitational back-reaction [43-47] a = eTGfi with e = 1, e = 10~^, 
and e = 10~^^ respectively. We have taken F = 50 everywhere. The dashed line shows 
the burst rate for loop sizes determined by the large scale dynamics of the network. We 
have taken a — 0.1 as suggested by recent numerical simulations [48]. Unfortunately, 
none of these models result in a signal detectable by LIGO. 
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Figure 6: Plot of JIqw as a function of the string tension for p — 10~^ and a — eTGji. 

The sohd hnes show the stochastic background generated at 75 Hz, corresponding to 
LIGO frequencies, the dotted hne shows the spectrum at 10^^ Hz, the band at which 
pulsar timing experiments are most sensitive. Again, unfortunately, none of these 
models result in a detectable signal. Although not shown here, large loops {a = 0.1) 
also do not result in a detectable signal. 
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4 Discussion and Conclusion 



We have calculated the gravitational waveform due to cosmic superstring reconnection, 
finding that bursts do not occur frequently enough to be seen with upcoming gravita- 
tional wave detectors, nor is the stochastic background strong enough to be detected. 
This conclusion remains true even for very auspicious values of currently unknown pa- 
rameters such as reconnection probabihty and loop sizes. Thus the most likely source 
of cosmic string gravitational bursts remains that resulting from cusps or kinks. 

This is doubly disappointing since one very attractive feature of the spectrum is 
that the gravitational radiation is polarized along the directions of {v, 6) at reconnection 
and so there must be a consistency relation with the probability distribution P{v,9) 
encoded in the kinks on horizon-sized strings. LISA, however, may be sensitive enough 
to detect bursts from reconnections and a detailed calculation including confusion noise 
from the network is underway. 

Despite this disappointing result from gravitational bursts, the techniques we have 
developed here can immediately be generalized to other (non-gravitational) particles, 
which could be observed in complementary experiments such as gamma-ray bursts 
or other astrophysical phenomena. Since the suppression is only power-law and not 
exponential in energy, there remains the possibility of producing very high energy 
particles. 
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A Appendix: Gravitational radiation formalism 

In the weak field limit the metric is assumed to be that of a perturbed Minkowski space 
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with \h^^\ <^ 1. In the harmonic gauge, defined by 

d,K - ld^h>; = 0, 

the hnearised Einstein equations read 

dxd%, = -167rG{T^, - ^r^^.T;) (A.21) 
where T^i, is the stress-energy tensor of the source. The retarded potential solution to 



Eq. (A.21) is given by 

r d^^' 1 
V(x,t) = 4G J ^— ^(r^,(x',t-|x-x'|)--r/^,r;(x',t-|x-x'|)). 

For sources localised in space and time we can construct the stress-energy tensor from 
the Fourier integral 

/oo 
duj T^,(x,cu)e-'"*, 
-oo 

with the inverse transformation given by 

1 r°° 

T^,(x,a;) = — / dt T^,(x,t)e*"*. (A.22) 
The retarded potential solution is 

dio / -—-iT,,i^',Lu) - -r/^,T;(x',^))e-^-(*-l^-^'l). (A.23) 

-oo J \ ^ 

We replace |x — x'| with r = |x| in the denominator of the integrand in Eq. (A.23) 

and in the exponent approximate |x — x'| ^ r — (1 ■ x'. This yields 

An poo /"I - 

/i^^x, t)-— J du e-*-(*-'-) J d^'x' (T^,(x', u) - -r7^.T;(x', a;))e-^'^"^-'''. 

If we now substitute for T^;^(x, c<j) using Eq.(A.22) we get 

^ J —oo 

"^^l^ dtj rfV (T,.(x',t) - ^r/,.T;(x',t))e^-(*-^-'). 
We can thus write the solution as the integral over the plane waves 



/oo 
du e^^{x,uj)e-'''-'' 
-oo 
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where 



1 



^ V . 2 



is the polarisation tensor, and 



oo 



T^^(A;) = — / dt I £x T^^(x)e^'=-^ 

271 J 

is the Fourier transform of the stress-energy tensor. 
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